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SUI€UARY 


A theoretical investigation was undertaken to 
evaluate the effect of small deviations from flatness 
on the behavior of simply supported square" plates under 
compression, it being aasvaned that the proportional limit 
is nob e3^oeded. This study involved the solution of the 
von Karm^ large-deflection equations according to the 
method used by Samuel Levy in NACA TN No, 81+6 . As a 
result of the investigation, it was fotnid that; 

(1) The effects of initial deviations from flatness 
upon buckle growth and effective width of plate are most 
marked at stresses around the theoretical flat-plate 
critical stress; at stresses well above or below the 
critical stress, the behavior of a plate with an Initial 
deviation from flatness is very much the same as that of 
an initially perfectly flat plate. 

(2) In two commonly used laboratory methods giving 
experimental critical stresses related to the start of 
rapidly increasing lateral deflections it may be expected 
that the larger the initial deviation from flatness the 
smaller is the experimental critical stress. 

( 3 ) The Southwell plot method of predicting theo- 
retical critical stresses for perfect specimens from 
experimental observations on actual specimens may not be 
expected to give, in general, satisfactory results when 
applied to flat plates, 
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lflTRODITCa?rON 


Plat and curved plates constl tut-e~basic structural 
elements in an all-metal airplane structure. Thie behavior 
and strength of both types of plate for various loading 
conditions have been treated in the literature in varying 
degrees . 

Ihe flat plate xinder compression has been treated 
very extensively in previous vork, but almost always with 
the restriction that the plate be perfectly flat. In 
ex'perimental studies of the compressive strength of stif- 
feners and stiffened panels that consist of flat-plate 
elements, the so-called flat platen have been found to 
behave not as perfectly flat-plates, for which buckles 
begin to form at a definite load called the critical load, 
but to behave rather as plates with slight deviations 
from, flatness, for which buckles begin tn grow with the 
beginning of loading. The rate at which the buckles 
grow with Increase in load is very slow at first but 
increases appreciably as the critical load is approached, 
after which the rate of buokle growth with, load gradually 
diminishes. 

Since every flat plate in practical use behaves as 
a plate with alight deviations from flatness, a theo- 
retical study of simply supported. square plates under 
compression was undertaken to evaluate the effect of 
small deviations from flatness on 

(1) Growth of buckles with increase in load 

(2) effective v/idth of plate 

( 3 ) Experimentally determined buckling or critical 

stress 

( 4 ) Applicability of the. Southv/ell plot method of 

predicting theoretical critical stresses for 
perfeotrly flat plates from experimental 
observations on actual plates 


The present study involved the solution of the von Karman 
large-deflection equations according to the method of 
Samuel Levy (reference 1) in order to obtain simulated 
test data for plates with several initial deviations from 
flatness. 
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The mathematical development of the equations required 
for the present investigation is given in appendix A for 
a rectangular^ plate under lateral pressure and longi- 
tudinal and transverse compression. In appendix B the 
special case of a square oiate subjected to compression 
in only one direction is treated. Throughout the investi- 
gation the assumption is made that the stresses remain 
below the proportional limit. 


SY?1B0LS 



D 


The coordinate system is shown in figure 10. 

plate length in x-direction 

plate length in y-direction 

effective width for load-carrying capacity 

effective width for stiffness against further 
edge compression 


flexural rigidity 


Et5 




E 

P 

k 


Young’s modulus 
stress function 

■positive integer used as subscript 
dimensionless deflection coefficient 



K. 


-’mn 


dimensionless coefficient for initial deflection 


Wr 


■^mn 


m,n,p positive integers used as subscripts 
Pj,g coefficients in infinite series for p^. 


normal ores sure 
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k 


q,r,s,t positive integers used as subscripts 


t 

U, VjW 

Wo 

Wmn 

Wq 

^>7>z 


xy 

Y’ 


2cy 


6 

6, 


e. 


S' » €\ 

X 


■or 


X 


a 

a. 


a 


plate thickness 

deflections of a point in the x-, y- , and 
z-dlrectlons , respectively 

initial deflection in z-directlon 

coefficients in infinite series for w 

coefficients in inflnl t-e~aeries for 

coordinates of a point measured from origin at 
corner of elate (see fig. 10) 

shearing strain 

shearing strain in middle surface 

total deflection at center of plate 

initial deflection at center of plate 

unit plate shortening in ^c-direotion (total 
shortening divided by a) 

unit plate shortening in y-direction (total 
shortening divided. by b) 

compressive strain in x- and y-directions , 
respectively 

middle -surface compressive strain in x- and 
y-directions , resoectively 

extreme-fiber compressive ! strain in x-direction 
at concave side of plate center 

extreme-fiber compressive strain in x-dlrectlon 
at convex side of plate center 

critical strain in x-direction for a flat plate 
subjected to compressive stress in x-dlrection 

Poisson’s ratio 

edge compressive stress in x-dlrection 
edge compressive stress in y-direction 
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average edge compressive stress in x-direction 

average edge compressive stress in y-directlon 

a compressive stresses in x-and y-directions , 

^ respectively 


p’x^ a’-n- middle-surface compressive stresses in x- and 
“ y-directlons , respectively 

extreme-fiber compressive bending stresses in 
X- and y-directions , respectively 



T 

T 

T 


xy 

t 

xy 

ti 

xy 


critical edge compressive stress in x-direction 


for a flat nlate 


/ 

V 


EiT^t^ ■ 


5(1-H^)b^ 


for a square plate) 


shearing stress 

shearing stress in middle surface 

extreme-fiber shearing stress due to bending only 


Deviations from Flatness Considered 

The types of initial curvature, deviations from flat- 
ness, selected for the numerical analysis are those com- 
ponents of the buckle pattern that first become predominant 
after the assumed-flat square plate buckles. 

Let the deflection surface of the simply supported 
square plate be represented by the equation 


m=l n=l 
m=l n=l 

where the absolute value of is the amplitude of a 
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component deflection surface 'witli m half waves In the 
x-direction and n half waves in the y-directlon, and 
is a nondimens ional measure of w^„ in terms of 

the plate thickness ~ * 

At zero load .the deflection surface w^ , which 
describes the Initial deviation from flatness, is given 
as 


w 


o 


CO CO 

t y~ y~ sin sin BEZ 

Z — °mn ^ ^ 

m=l n=l 


( 


Kr 


•'mn 


0 for perfectly flat plate 


) 


Ylian a simply supported, oerfectiy fLat square plate 
first buckles, the deflection surface w is of the form 


tK,, sin ain BZ 
b b 

but as the deflection increases, other components appear, 
the most Important of which is 


tK,, am iis sin 22 
3^ b b 

It is therefore reasonable to assume that 
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sin 2^ sin, 2:1 
b b 
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and 


tK„ aln sin ^ 
fc D 


are the two moat i mportant t .To ea of in itial deviation from 
flatness;' eTnce the first of these types has^a'- greater 
errect. r.nan the second, more variations were studied of 
the first than of the second. The particular combinations 
of Initial curvature of the ^02^2 ^°31 con- 

sidered -for numerical computation are given in the fol- 
lowl tig~tra bl'e l ^ ^ 


°11 


0 

0 (flat) 

.01 

0 

• 0l4. 

0 

.10 

0 

0 

.01 ■' 

.OI4. 

.01 , 


RESULTS AND DISCUSSION 


The assumption is made in the analysis that all edges 
of the plate remain straight -vdaen buckling occurs, and 
that tlie side edges are allowed to translate in the plane 
of the plate to such a position that the resultant of the 
transverse forces on the side edges of the plate is zero. 
The edge positions and stress distribution corresponding 
to such an assumption for an initially perfectly flat 
date are Indicated in figiire 1 » vhen the plate Is not 
perfectly flat the stress distributions are, from the 
beginning of loading, more or less of the general type 
Illustrated by the "after- buckling" case shown in figure 1 . 

The theoretical study of the effect of deviations 
from flatness on 


( 1 ) Growth of buckles with increase in load 
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(2) 'Effective width of plate 

( 5 ) Experimentally defined buckling or critical 
stress 

(if.) Applicability of Southwell plot method of- 

predioting theoretical critical stresses for 
perfectly flat plates -from exoerlmental 
observations on actual plates 


disclosed relatlonshins which would be of most general 
value when presented in nondimenslonal fom by the use 

t ’ o e b 


of such natural ratios as 


b* 


and 


e 


Poisson’s’ ratio p, is 


taken as 


X 

0.516 


for the 


comnutations . The comcuted results apply only if the 
proportional limit is not exceeded. 


Growth of Buckles with Increase in Load 

The way in which the net center deflection 6-6^ 
JLnoreases with load as measured by the average stress 
is shown in figure 2. ’'hen a perfectly flat-plc.te 
is sub,1eoted to an increasing edge compression, no 
deflection of the date out of its original plane occurs 
until a definite stress, called the critical stress, is 
reached. An increase in load above the critical load 
causes large deflections which, grow less and less 
raoidly as the load increases more and more. (See 
curve for K_ = = 0 in fig. 2 . ) 

"’ll °^1 ; 


’^toen a d^te with an Initial de\d.ation from flatness 
is subjected to compression, the deflections that result 
frcm application of load grow very slowly with the first 
increments of load but grow very rapidly as the critical 
load is approached; near the critical load the rat^ of 
increase of 6 - 6g with further increase In load 
diminishes and 6 - 6 q aoproeohes the deflection for the 
perfectly flat plate at values not far above the critical 
load, {See curves other than =0 in fig* 2.) 


’ll 


°51 
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Effective width, of Plate 

Ihe f\mdamental information for evaluating the 
effect of deviations from flatness on two t;ypes of effec- 
tive width of date is presented in figure 
ratio of average stress to critical stress is 

plotted against the ratio of xmlt plate shortening to 
critical strain ^a/^cr • 

The usual "effective width" bg, which is associated 
with the load-carry! ng _. canac.it y of the plate, may be 
defined (see last paragraph of appendix B for derivation) 
by the equation 


b 


e 


'a/'or. 


Thus the ratio of the effective width for load-carrying 
capacity to the actual width of . the plate is equal 

to the ratio of ordinate to abscissa in figure 5* 
figure Jx the effective-width ratio tg/b for load- 
carrying capacity is plotted against the ratio of average 
stress to critical stress for plates with dif- 

ferent amounts of initial deviation from flatness. 

Another tyoe of "effective width" b’ is associated 
v.d hh the stiffness gf ^^^ the plate a gainst fur the jc- compr e-s— 
sTon ^that is, dd^yd^gj^ and is defined by the equation 


(A. detailed discussion of this type of effective width 
is contained in reference 2 where it is referred to as 
"reduced effective width.") Thus the ratio of effective 
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width for stiffness to actual wi.dth b'e/^ given 
directly by the slope of a curve of ^a/'^cr against 

^aAcr • Such curves are shown in; f 3 gure 5 for plates 

with different" amoohts of initial deviation from flat- 
ness. In figure 5^ effective-width ratio for stiff- 
ness is plotted against the ratio of average 

stress to critical stress . 

O X 

'die effects of initial deviations from flatness upon 
the effective width of simply supported square plates arc 
most marVced at stresses around the theoretical flat-plate 
critical stress. The effective width for load-carrying 
capacity of a olate with an initial deviation from flat- 
ness is, at all values of stress, less than that of an 
initially perfectly flat plate (fig. 4). The effective 
width for stiffness against further compression, however, 
is less than that of a flat plate at stresses below the 
theoretical critical stress but greater at stresses above 
the critical (fig. 5)* At stresses well above or below 
the critical stress the effect of initial deviations upon 
either effective width of plate becomes negligible. 


Experimentally Defined Buckling 
or Critical Stress 

Since buckles begin to grow with the beginning of 
loading for a plate that has an initial deviation from 
flatness, there can be no buckling stress for an actual 
olate in the strict theoretical sense; hov/ever, just as 
a defined yield stress has been found useful for 
materials that have no actual yield stress, so a defined 
buckling or critical stress for a plate can convey much 
meaning to a practicing engineer. 

A perfectly flat plate" remains flat until the 
buckling stress is reached, whereupon buckles are produced 
that grow rapidly with very small_further increases in 
load. (See curve labeled = K_ =0 in fig. 2.) 

Oil °51 

Since the rapid increase of lateral deflection with load 
is an important aspect of the buckling of a perfectly flat 
plate, it Is reasonable to retain this concent as the 
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basis of an experimentally determined buckling stress for 
an actual plate. Because, however, strains are usually 
more easily measured than deflections with the equipment 
available at present and because the manner of growth of 
the extreme-fiber strains at either side of the buckle 
crest is indicative of the increase In lateral deflections, 
these strains instead of the lateral deflections are 
frequently m.easured in determining the buckling stresses 
of plates. 

In two NACA papers on the compressive strength of 
stiffened panels (references 5 s- method that will 

be referred to herein as the "strain-reversal method" has 
been used to obtain experimental buckling stresses. Ihe 
critical stress obtained by this method is defined as 
the stress at which the extreme-fiber strain 

convex side of the buckle crest stops increasing and begins 
to decrease. Critical stresses obtained according to 
this definition are indicated bj circles in figure 6 , 
where the calculated extreme-fiber strains in the direc- 
tion of loading at either aide of the buckle crest are 
plotted as the abscissa against the average applied stress 
as the ordinate for plates with various assumed initial 
deviations from flatness. 

In a series of NACA papers on the plate compressive 
strength of metals (references 5 10 ), the experimental 

buckling stress was determined by a method that will be 
referred to herein as the "top-of-the-knee method." The 
critical stress, according to this method, is essentially 
the stress corresponding to the top of the knee of a 
curve of stress against lateral deflection. If the lateral 
deflections cannot be readily measured, any other quantity 
that increases in substantially the same manner as the 
lateral deflections may be plotted instead. One such 
quantity is difference in strains in the 

direction of loading at the two sides of the buckle crest. 
Buckling stresses obtained by the top-of-the-knee method 
are indicated by squares on the curves of figure 7 for 
the several assumed Initial deviations from flatness 
considered in the present paper. 

The critical stresses obtained by the two methods 
,iust described are shown transferred to the stress- 
deflection curves of figure 8 . Both methods give critical 
stresses that lie roughly on the knee of the stress- 
deflection curve, the strain-reversal critical stress 
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being the more conservative of the two. Since the knee 
of the curve may be jtidged to be the transition stage 
from a low to a high rate of increase of lateral deflec- 
tion wJ.th load, both methods can be considered appro- 
priate- for obtaining experimental buckling stresses, 
according to the criterion adopted herein as the basis 
for defining an experimental buckling stress. Regardless 
of which method is used, a certain degree of personal 
judgment- must be exercised in the selection of the buckling 
stresses, and the experimentally obtained buckling stress 
becomes more and more uncertain as the initial deviation 
from flatness increases. Despite this uncertainty, it 
is clear that the effect— of initial deviations is to 
reduce the experimentally determined buckling stress. 


Critical Stress and Postbuckling 
Behavior of Plate 

The experimental definitions of critical stress just 
discussed were based on the concept of rate of growth of 
lateral deflection with load as the theoretical critical 
stress is aporoached. Since the rate of growth of lateral 
deflections below the theoretical critical str-ess is 
strongly influenced by the initial deviation from flatness, 
it is natural and perhaps desirable that any method for 
determining experimental critical stresses should be 
sensitive to initial deviations from flatness. 

At stresses considerably beyond the buckling stress, 
however, the behavior of the plate is essentially inde- 
oendentof initial deviations from flatness. (See 
figs. 2 to 7*) describing the. elastic plate behavior 

in this range of stress, it is therefore advantageous to 
use, if possible, the theoretical crj.tical stress, which 
is also independent of the initial imperfections of the 
plate. 


Southwell Plot Method of Predicting 

Theoretical Critical Stress 

Because the critical stresses for perfect specimens 
are valuable both for checking theories and predicting 
oostbuckling behavior and because in many cases such 
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critical stresses are difficult to compute, a method of 
determining them from experimental data on imperfect 
specimens appears to he very desirable. Southwell has 
produced such a method for columns (reference 11), which 
has been modified in reference 12 for easier application 
to experim.ental work. Although this method applies to 
flat plates under compression according to smail-def lection 
theory (o. 521 of reference I 5 ) it has been found some- 
what unsatisfactory for general use in the laboratory. 

An essential requirement for the success of the 
Southwell method is that at least a part of the stress- 
deflection curve approximate a rectangular hyperbola 



(for simply supported plates) 



(for pln-endod columns) 


where 6 is the total deflection of the center of the 
col\min or plate, 5^ is the initial deflection of ^he ' 
center of the column or plate, a is the average compres- 
sive stress on a column, and O^p is the critical com- 
pressive stress. The test data from which the Southwell 
plot is constructed must then be taken from this part of 
the curve. Tn the light of this requirement, the reason 
for the success of the Southwell method when applied to 
columns and for its failure when applied to plates is 
made clear in figure 9 j where the stres-s-def lection 
curves for simply suoported columns and plates are com- 
pared with the rectangular hyperbolas that they must 
approximate. (The stress-deflection curve lii fig. 9(^) 
is for a column of vjhich the initial shaoe is a circular 
arc and was calculated by large-deflection theory by a 
method similar to that used on pp. 69 to 72 of reference I 5 .) 
Tn the case of a column, regardless of the initial shape, 
the requirement that a part of the stress-deflection curve 
approximate a rectangular hyperbola is certain to be met 
as the stress approaches the critical value. In the case 
of a plate, however, the stress-deflection curve usually 
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does not ao'Droyiinate a reotangulsr hyperbola at stresses 
near the critical stress unless the initial deviation 
from flatness haooens to be very small (fig. 9(s-))» 
Exoerience in the structures research laboratory of the 
NACA Langley Laboratory indicates that the initial devia- 
tion from flatness and the inevitable inaccuracy of 
loading are usually not sufficiently small to guarantee 
this close approximation. Even if these factors are 
sufficiently small, difficulty is encountered in measuring 
with sufficient accuracy the very small strains or deflec- 
tions used to make the Southwell plot for this special 
case. 


The close agreement between the plate s tress - 
deflection curves and the rectangular hyoerbolas at low 
values of stress (fig. is largely fortuitous and 

occurs because the type of initial deviation from flat- 
ness assumed is the same as the lowest buckling mode. 


OONCLtrSTOKB 


The following conclusions are based on an analysis 
of simulated test data calculated by large-deflection 
theory and apoly to the elastic behavior of simply 
suoported square plates under compression. 

1. The effects of initial deviations from flatness 
upon the bucl<le growth and effeotive width of simply 
supported square plates are most marked at stresses 
around the theoretical flat-olate critical stress. At 
stresses well above or below the critical stress the 
behavior of a elate with an Initial: deviation from flat- 
ness is very much the same as that of an initially 
perfectly flat plate. 

2. Tlie effective width for load-carrying capacity 
of a plate with an initial deviation from flatness is, 
at all values of stress, less then that of an initially 
perfect plate. The effective width for stiffness against 
further compression, however, js less than that of~ a flat 
plate at stresses below the theoretical critical stress 
but greater at stresses above. 

5 . Either of two commonly used laboratory methods, 
the strain-reversal method and the top-of- the-knee method, 
gives experimental critical stresses that are lower than 
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the theoretical ' flat-plate critical stress, the strain- 
reversal trethod generally giving more conservative values 
than the top-of-the-\fnee method. Such experimental 
critical stresses are significant as sin indication of the 
start of rapidly growing lateral deflections. The effect 
of increasing the initial deviation from flatness is to 
reduce the experimental critical stress as determined by 
either method. 

k. The Southwell plot method of predicting theoretical 
critical stresses for perfect specimens from experimental 
observations on actual specimens may not be expected to 
give, in general, satisfactory results when applied to 
flat plates . 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. , May I 5 , 19^6 
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APPENDIX A 

SOLUTION OP VON KARKAN Ei^UATIONS 


I'he differential equatlona derived by von ^finiiAn for 
the equilibrium of a flat rectangular plate under the 
action of normal pressure p„ and central surface stresses 
described by the stress function F are as follows (refer- 
ence 1) ; 


dx^dy^ dy^ 

d^w . ^ d^w , d^w 
dx ' dx^dy^ d y^ 


d^w d^w 
dx^ dy^ 


Pz 

D 



(Al) 


^d^F dfw dfp pdfp. 

vdy^ dx^ 


dx^ 67^ 



(A2) 


(The coordinate system is shown in fig. 10.) 

If the plate has an initially curved surface 
described by the function w„, the strains used to 


establish equation (Al) must 


Strains 

Flat plate 

-e ^ 

^ dx 

^ idjf 


dy 



-V ^ 

dv dw 

♦ — - - 1 ■ — 

dw 

^xy 

d X dx 

dy 


be revised as follows: 
Initially curved plate 

^ 

dx 2^x/ 2ydx / 

^ + U^f . i/^Y 

dy 2^jy ~ ^\b7 ) 
du ^ ^ ^ d Wq dWp 

dy dx dx dy ~ dx 6y 
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The left side of equation (A2) Is obtained from 
expressions for bending and twisting moments; and, since 
the moments depend not on the total curvature but only 
on the change In curvature of the plate, the net deflec- 
tion w - Wq Instead of the total deflection w should 
be used on the left side for an Initially curved plate. 
The right side of equation (A2) gives the combined effect 
of the lateral load P 2 and forces In the plane of the 
plates; and, since the effect of the forces In the plane 
of the plates depends on the total deflection w, the 
right aide of equation (A2) remains unchanged for an 
initially curved plate. 


The equations of equilibrium for an initially cxirved 
plate therefore become 


^ + 2 ^ = -p 


\hx6yy 




62 


6^Wq 62wq 




(A5) 


- Wq) 6^(w - Wq) 5^(w - Wq) 

— + 2 r 1 + 


6x^ 


6x2by2 


6y^ 


Ez 

D 


i. { 6^w ^ 6^F 6^w ^ 62p 

^ ^ \6y2 6x2 ^ ^^2 ^y2 6x6y 6x6y> 


( Aij.) 


The middle -surface stresses are 

q2f 




0^=-^ 
^ 6x2 

- 6^F 

‘ 6x 6 y 


> 


(a5) 
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The extreme -fiber bending and, shearing stresses due 
to bending alone are 


„ ^ Et - Wq) 6^{w - Wq) 

c;" = + ; — + IJ, 

^ 2(1 - IJ,^) dx^ 


■ " = + 


Et 6‘^(w - Wq) ^ 6*^(w - Wq) 

2 ( 1 - ) 6 6 x^ 


(A?) 


- + Et - Wq) 

2(1 + [j.) dxdy 


A solution of equations (A5) and (Aii) must satisfy 
the boundary conditions that the deflection and the 
bending moment per unit length along the edges be 
zero; 
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At X = 0 and x = a, 


- Wq) - Wq) 

^ _Q 


6x" 


6y2 


At 7 = 0 and y = b. 


6^(w - Wq) 5^(w - Wq) 




+ \i- 


6x2 


-^ = 0 


These conditions are satisfied by the two Fourier 
series 


w 


Z Z ^ 




00 05 

= ZZ 

m=l n=l 


sin ^ sin Sm 
°!nn a b 


(a8) 


The normal pressure may he expressed by the Fourier 
series 


05 C5 

='^ sin SE 3ln 2|2 (AQ) 

r=l s=l 


By substitution, equation (AJ) is found to be 
satisfied if 


p = - 


>y^ CJx,x2 


p=0 q=0 


' ■ ‘’Ppq) 


prrx 


£2IZ 


cos cos (AlO) 

a b 
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where and are the average applied stresses In 

the X- and. y-dlrectlons , respectively, and where 


Q 



and 

[kt(p - k)(q - t) - (q - t)^] '^(p-k)(q-t) 

k=l t=l 

if q ^ 0 and p. ^ 0 

B-j_=0, if p = 0 or q=:0 
CO q~i 

Bg = ^ ^ [kt(k + p)(q - t) + (q - t)^] W(k+p)(q_t) 

k=l t*l 



If q 0 




"3 


= 0, if q = 0 

00 q-l 

[(k + p) tk (q - t) -H (k + p)^ (q - t)^ 

k=i t=l 



^\k+p)t \(q-t) 


if q 7^ 0 and p / 0 


B^=0, if q = 0 or p=0 
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P-^ 


= \ \ [kt(p - k)(t + q) + (t + q)^] Wj^.^ v<(-p_ 


k=l t=l 


k) ( t+q) 


if p 7 ^ 0 


= 0 , if p = 0 ■ 
p~l ” _ 

B5 =y~ \ L^(t + q)(p - k)t + k^t^J "^(p-k)t 

k=l t=l 

if q 7 ^ 0 and p 7 ^ 0 


0 


= 0 , if q = 0 or p = 0 


OP OP ^ 

B6 = + q) - + q)^J 

k=l t'5'1 


if q 7 ^ 0 or p 0 


3 ^ = 0 , if q = 0 and p = 0 


CO 00 

B? "" ~ ""(k+p)t 

k=l t=l 

if q 7 ^ 0 and p 7 ^ 0 
= Oj if q = 0 or p = 0 
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a> (■"> 



k=l t=l 


[(k + p) tk(t +-g) - (k + p)^ (t + q)^l 


^{k+'o)t ^k( t+q) 


Bs = 0, If 


00 OO 



k=l t=l 


3^ =0, if 


if q 0 and p ^ 0 
p = 0 or q = 0 

[(k + p)(t + q)kt - (k + p)2 t^] 

If V ^ 0 or q 7 ^ 0 
p = 0 and q = 0 ; 

q 



= E ^ 

11(9^1 + q^f) 



The terms are the same as the B terms except 

that w terms are replaced by Wq terms. 

The terms bQQ and ^oqq iJ? ay have any arbitrary 

value but they are -of no importance because they vanish 
In the oartial derivatives of the stress function F 
aouearing in equations (A3) and (Aif.). 
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Equation (aU) is satisfied if 


vvhere 


rs 


= D(\ 


'N /r^2_ + 


rs "Ox.s;\^ a2 


- Og^tw^^gr^ h| - OT^tw^g! 
BT 


trA 


(^rs ~ ^Ops) 


^rs ^ ^ 

n=i 


9 

Aq ~ Aq 

"^rs / '^n 

n=l 


Ai = - 


r s 2 

y~ "y [(s - t)k - (r - k)t] 

fel t=l 


03 CO 


A^ - - 



[t(k + r) - k(t + s)3 b 


kt 



(All) 


-k) (s-t) ^kt 


^(k+r) ( t+s) 


k=0 t=0 
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The e quat?! o ns f o r 
fo3? 5 


A^ . . . A^ are the same as those 
respectively, except that the h 


terras are replaced by b^ terms. 


By integration the function P is found to satisfy 
the conditions 





b 


for X = 0 and x = a and 



for y = 0 and y = b 

The remaining conditions thet the solution must 
satisfy are that the edges x=0, x=a, y=0, 
and y = b remain straight. These conditions require 
that a£g^ be independent of y and that be^ be inde- 
oendent of x. By actual evaluation, it is found that 



These expressions for and are independent of 

both X and y. 
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APPENDIX B • 

APPLICATION TO SQUARE PLATE 


The formulas required for numerical v,’ork are given 
for the special case of a square plate subjected to 
compression in only one direction. Let 


a = p 



p, = 0,516 


The system of equations (All) becomes simplified in this 
case to 



When symmetrical initial curvature and syraraetrical 
buckling are assumed, eind when only the first six com- 
ponents of the double harmonic series with odd subscripts 
are considered, may 

be deterralned from six equations of the type of equar 
tion (All) when ^a/^cr initial curvature are 

known. These equations are similar , to those in table II 
of reference 1. 
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Extreme-fiber strains .- From equations (a 6), (A7)# 
and (AlO), the extreme-fiber strains at the plate center 
are given by 




+ 3(1 - 


p=0 q=0 


Et 




2+3. 

2 




00 CO 


- 3ti(l - 



p*0 q=0 


2 211 
(V ■ SO*'"’ " 


OO 00 


+ 1,5(1 - n^) 


(-1) = 


m=l,3,5 n=l,5,5 


>(B2) 


^2 


CO CO 




+ 3(1 - |J-^) 


p«Q q=o 


-i 

Et 


2 


Eta 

2 


CO CO 


-5fi(l-p2) 


Et2 


(V - »0p3)<-« 


Eii 

2 


p=0 q=0 


CO CO 


- 1.5(1 -4x2) 




in=l»3»5 51=1,3,5 


The difference of the extreme-fiber strains 
in the x-directlon at the center of the olate divided 
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by e Is given by the followirjg equation; 

Lr X 



00 



m=*l,5,5 


CD 








Lateral deflection .- The deflection at the oenbor 
of the olate divided by the plate thickness t is 


OO CO 



m=l,3,5 n=l,3,5 

The net deflection at the 
by the plate thickness t is 


CO 00 



m-'L,3,5 n=l,5,5 



center of- the plate divided 



Effective width .- If a flat plate shortens v/ithout 
buckling, the average stress is EF_. The average stress 

Bi 

on a buckled plate or a pl^te vd th initial curvature is 

which is less than SCg. In other words, the buckled 

plate or the plate with an initial deviation from flat- 
ness is equivalent in lo ad -carr-ying capacity to a flat 
and unbuckled plate of the same thickness but smaller 
width, '.{he width of this hyootheticaL plate, or the 
effective width b^ of_the actual plate, can be obtained 
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by equating the load-carrying capacities of the two 
dates. -Thus 




Since, for elastic buckling. 


hence 


E 


^cr 


cr 


X 

X 



(b6) 
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Fig. la-c 



(a) Before loading. 


(b) After loading; just 
before buckling. 


(c) After buckling. 


NATIONAL ADVISORY 
COMMITTEE FOR AERONAUTICS 


Figure 1«- Edge-stress distribution for a simply 
supported flat plate under edge compression. 



Fig. 2 
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Figure 2,- Variation of net-center-deflection ratio 

with average edge-compressive-stress ratio for simply 
supported square plates with slight "initial deviations 
from flatness. 
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Scrx 

Figure 5,- Variation of unit-plate-shortenlng ratio 

with average edge-compressive-stress ratio for simply 
supported square plates with slight initial deviations 
from flatness. 
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Fig. 5 
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Figure Variation of effective-width ^ratio for stiffness 

with average edge-compresslve-stress ratio for simply 
supported square plates with slight initial deviations 
from flatness. 



Fig. 6 
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Figure 6.- Variation of ejctreme-f iber-strain ratios at 
plate center with average edge-compresslve-atress ratio 
for simply supported square plates with slight initial 
deviations from flatness. 









Fig. 8 
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Figure 8.- Comparison of buckling-stress ratios obtained by 
strain-reversal method and top- of- the- knee method. ^ 
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Fig. 9 
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(b) Rectangular column 
with initial shape of 
circular arc. 

Length = JOt. 


s 
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Figure 9»-^ Comparison of calculated stress-deflection 

curves for simply supported columns and plates with the 
rectangular hyperbolas that the curves must closely 
approximate if the Southwell plot method is to apply. 




Fig. 10 
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Figure 10.- Coordinate system ibr rectangular plate. 



